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We study the influence of the regular component of the Galactic magnetic field (GMF) on the 
arrival directions of ultra-high energy cosmic rays (UHECRs). We find that, if the angular resolution 
of current experiments has to be fully exploited, deflections in the GMF cannot be neglected even 
for E = 10^° eV protons, especially for trajectories along the Galactic plane or crossing the Galactic 
center region. On the other hand, the GMF could be used as a spectrograph to discriminate 
among different source models and/or primaries of UHECRs, if its structure would be known with 
sufficient precision. We compare several GMF models introduced in the literature and discuss for the 
example of the AGASA data set how the significance of small-scale clustering or correlations with 
given astrophysical sources are affected by the GMF. We point out that the non-uniform exposure to 
the extragalactic sky induced by the GMF should be taken into account estimating the significance 
of potential (auto-) correlation signals. 



PACS numbers; 98.70.Sa, 98.35.Eg 

I. INTRODUCTION 

Despite of more than 40 years of research, the field 
of ultra-high energy cosmic ray (UHECR) physics poses 
still many unsolved problems [l| . One of the most impor- 
tant open issues is the question at which energy astron- 
omy with charged particles becomes possible. The an- 
swer depends both on the chemical composition of UHE- 
CRs and on the strength of the Galactic and extragalac- 
tic magnetic fields. Extensive air shower experiments 
can in principle measure the chemical composition of the 
CR flux. However, the indirect measurement methods 
and the differences in the predictions of hadronic inter- 
action models make the differentiation between proton 
and heavy nuclei primaries a theoretical and experimen- 
tal challenge 0. Other signs for proton or nuclei pri- 
maries are therefore highly desirable. 

Complementary information on the charge of the pri- 
mary may be obtained by studies of the arrival directions 
of CRs. Historically, effects of the geomagnetic field on 
the CR propagation were decisive to understand the na- 
ture of low-energy CRs: The discovery of the latitude 
effect proved that a significant fraction of cosmic rays 
is charged, and the east-west asymmetry demonstrated 
the predominance of positively charged primaries It 
is natural to ask if the weaker magnetic fields known 
to exist on larger scales like the Galactic magnetic field 
(GMF) might play a similar role at even higher energies, 
thus providing important information about the charge 
composition and the sources of the UHECRs. 

A possible signature of proton primaries is the small- 
scale clustering of UHECR arrival directions. The small 
number of sources able to accelerate beyond 10^^ eV 
should result in small-scale clustering of arrival direc- 
tions of UHECRs, if deflections in magnetic fields can 
be neglected or, to some extent, even in presence of a 
magnetized medium (see e.g. Q). For nuclei with higher 
electric charge Ze, the deflections in the GMF alone di- 
lute a small-scale clustering signal even at the highest en- 
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ergies observed. Therefore, the conflrmation of the small- 
scale clustering observed by the AGASA experiment at 
energies above 4 x 10"'^^ eV would favor the hypoth- 
esis of light nuclei primaries, in particular protons. At 
present, the statistical signiflcancc ascribed to the clus- 
tering signal varies strongly in different analyses 0, |^ . 
Moreover, the HiRcs experiment Q has not conflrmcd 
clustering yet, but this finding is still compatible with 
expectations [3, W, ■ The preliminary data of the Auger 
Observatory have been searched only for single sources, 
with negative result [To| . 

Main aim of this work is to quantify the effect of the 
Galactic magnetic field (GMF) on the arrival direction of 
UHECRs and to study the possibility that corrections for 
the GMF role may provide important information about 
the charge composition and the sources of the UHECRs. 
In particular, the impact of the regular GMF on the clus- 
tering signal and on correlations of UHECRs with BL 
Lacs is addressed. Since UHECRs have presumably ex- 
tragalactic origin, the possibility of UHECR astronomy 
relies on the negligible effect of extragalactic magnetic 
fields. Existing simulations agree on several qual- 
itative features, but disagree on the magnitude of the 
UHECR deflections. In the following, we assume opti- 
mistically that such deflections are small in most of the 
extragalactic sky. In Sec.^ we review the main features 
of three GMF models presented previously in the litera- 
ture. In Sec. mil wc discuss the role of the GMF for the 
propagation of UHECRs, and the method we use to as- 
sess the significance of a possible small-scale clustering in 
UHECRs data. In Sec. II VI we apply these concepts to the 
AGASA data set of events with energy E >Ax 10^^ eV, 
first to autocorrelation and then to test suggested cor- 
relations with BL Lacs. In Sec. we summarize our 
results and conclude. 
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II. GALACTIC MAGNETIC FIELD MODELS 

The first evidence for a Galactic magnetic field was 
found more than 50 years ago from the observation of 
linear polarization of starlight Meanwhile, quite de- 
tailed information about the GMF has been extracted 
mainly from Faraday rotation measurements of extra- 
galactic sources or Galactic pulsars However, at 
present it is not yet possible to reconstruct the GMF 
solely from observations (for an attempt see 0| ) , and in- 
stead we will employ three phenomenological models for 
the GMF proposed in the literature. The GMF can be 
divided into a large-scale regular and a small-scale tur- 
bulent component, with rather different properties and 
probably also origin. The root-mean-square deflection 
<5rms of a CR traversing the distance L in a turbulent 
field with mean amplitude Srms is (e.g. [Tsj ') 
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where Lc denotes the coherence length of the field, £'20 
is the energy in units of 10^° eV, and L ^ Lc has been 
assumed. Recently, Tinyakov and Tkachcv noted that 
the latter condition could be not fulfilled, at least for 
some directions in the sky However, their analysis 
based directly on the observed turbulent power spectrum 
confirmed that the deflections in the random field are 
typically one order of magnitude smaller than those in 
the regular one. Therefore, we shall neglect the turbulent 
component of the GMF in the following. 

The regular GMF has a pattern resembling the one 
of the matter in the Galaxy and has different properties 
in the disk and the halo. In the disk, the field is es- 
sentially toroidal, i.e. only its radial {Br) and azimuthal 
{Bg) components are non- vanishing. The disk field can 
be classified according to its symmetry properties and 
sign reversals: antisymmetric and symmetric configura- 
tions with respect to the transformation of the azimuthal 
angle 9 ^ 9 + tt are called bisymmetrical (BSS) and ax- 
isymmetrical (ASS), respectively. According to the sym- 
metry property with respect to a reflection at the disk 
plane {z -~z), the notation A or S is used: in the first 
case, the field reverses sign at z = (odd field), while 
in the second case it does not (even field). Theoretical 
motivations and observations in external galaxies |l7l| as- 
sociate the presence of field reversals far away from the 
Galactic center (GC) to a BSS geometry: In our Galaxy, 
there are probably 3-5 reversals. The closest one is at a 
distance of 0.3-0.6 kpc towards the GC, where the higher 
values seem to be confirmed from the new wavelet data- 
analysis used in [l^ IT^. and about 0.6 kpc is the value 
suggested in the review Moreover, there is increas- 
ing evidence for positive z p arity (configuration S) of the 
GMF near the Sun [11 El 113 . 

In Galactocentric cylindrical coordinates, the field 
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FIG. 1: The Galactocentric frame used in this paper, together 
with the Solar position (circle) along the y-axis and the ori- 
entation of the Galactic polar angle 6. The corresponding 
Galactic longitudes are also shown, as well as the Galactic 
spiral arm model as given in [2H . 



components in the disk can be parameterized as 
Br ^ B{r,9)smp, Be = B{r,9) cosp , 



(2) 



where p is the pitch angle and Rq ~ 8.5 kpc is the Galac- 
tocentric distance of the Sun, cf. Fig. ^ Estimates for 
the pitch angle are p = —8° ± 2° from pulsar and 
starlight polarization data, but other observations point- 
ing to a value of p ~ —13° — 18° also exist [l9j . 

The function B[r,9) is traditionally modeled reminis- 
cent of the spiral structure of the matter distribution in 
the Galaxy as 



B(r, 61) = b{r) cos 61 - 
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In terms of the distance d to the closest sign reversal, ^0 
can be expressed as ^0 = {Rq + rf) cxp(— ^ tanp). The 
radial profile function h{r) is generally assumed to fall off 
as (X 113,01 1 consistent with pulsar measurements 
113 ■ The behavior of the disk field in the inner region 
of the Galaxy is less known, but clearly the field can not 
diverge for r — > 0. For r > J'max, tbc field is turned off. 
In the following, we will fix Tmax = 20 kpc. The vertical 
profile of the field outside the plane z = is modeled by 



B{r,9,z)^ f{z)B{r,i 



(4) 



Despite remaining uncertainties, the regular magnetic 
field in the thin disk is yet much better known than other 
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components, namely the halo (or thick disk) field and 
a possible dipole field. The first one could dominate at 
large Galactic latitudes and the second one may be of cru- 
cial importance near the center of the Galaxy. Because of 
the huge volume occupied by the halo field, it may play a 
determinant role for UHECR deflections, while the pos- 
sibly much higher strength of the field in the center of 
the Galaxy might prevent us to access some directions in 
the UHECR extragalactic sky (see Sec. 

For the halo field, an extrapolation of the thin disk 
field into the Galactic halo with a scale height of a few 
kpc has often been assumed (e.g. |23,I3)- This minimal 
choice is in agreement with radio surveys of the thick 
disk and mimics the expected behavior of a "Galac- 
tic wind" diffusing into the halo. However, Faraday rota- 
tion maps [2II23 of the inner Galaxy (-90° < Ig < 90°) 
and of high latitudes (|t»G| > 8°) favor a roughly toroidal 
component in the halo, of opposite sign above and below 
the plane (odd z parity or configuration A) and with an 
intensity of 1-2 ^G |23|. Moreover, there is some evi- 
dence for a Bz component of about 0.2 at the Sun 
distance |29| that could derive from a dipolar structure 
at the GC [301 ■ In the filaments already detected, the 
field strength almost reaches the mG scale Even 
if this intriguing picture is roughly consistent with the 
one expected if a AO dynamo mechanism operates in the 
Galactic halo, it needs observational confirmation. For 
example, there is no general consensus about the exis- 
tence of such high-intensity magnetic fields in the central 
region of the Galaxy, see e.g. [3lll33 |. 

In the following, we review three phcnomcnological 
models that parameterize the regular GMF. These mod- 
els are characterized by different symmetries, choices of 
the functions b{r) and f{z) and parameters. 



A. TT model 

Tinyakov and Tkachev (TT) examined in if cor- 
relations of UHECR arrival directions with BL Lacs 
improve after correcting for deflections in the GMF. 
They assumed b(r) oa for r > rmhi — 4 kpc, and 
b{r) = const, for r < r^in- The field b{r) was normalized 
to 1.4 fiG at the Solar position. The pitch angle was cho- 
sen as p — —8° and the parameter d fixed to —0.5 kpc. 
They compared a BSS-A and a BSS-S model and found 
that for the former model the correlations with BL Lacs 
increased. This model has an exponential suppression 
law. 



f{z) =sign(z)cxp(-|z|/zo): 



(5) 



with zq = 1.5 kpc chosen as a typical halo size. No dipole 
component was assumed. 



B. HMR model 

Harari, Mollerach and Roulet (HMR) used in js^l a 
BSS-S model with cosh profiles for both the disk and 
the halo field with scale heights of zi = 0.3 kpc and 
Z2 — 4: kpc, respectively, 



1 



+ 



2cosh(z/zi) 2cosh(2;/z2) ' 
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Thus the disk and halo field share the same spiral-like 
geometrical pattern. The function b{r) was chosen as 
b{r) = 3Ro /r tanh'^ {r / ri) fiG with ri = 2 kpc, hence re- 
ducing to 6(r) oc r^^ for r 3> while vanishing at the 
Galactic center. The pitch angle was fixed to p = —10°, 
and ^0 = 10.55 kpc. This model represents a slightly 
modified and "smoothed" version of the BSS model dis- 
cussed by Stanev in Apart for the vertical profile 
/(z), the main differences with respect to the TT model 
are the z parity and the r — > behavior of the field. 



C. PS model 

In Prouza and Smida (PS) used for the disk 

field the same BSS-S configuration as with a sin- 
gle exponential scale height zq and b{r) as described in 
Sec. Ill AI In the slightly modified version we use here, we 
fix zo = 0.2 kpc, p = —8°, d = —0.5 kpc and normalize 
the local field to 2 iiG. Apart for the larger field-strength, 
the main difference with the TT model is the parity of 
the disk field, which we take here to be even as in |34l |. 
Additionally, we consider a toroidal thick disk/halo con- 
tribution. 



Btx = —Bt sign(z) cos 9 , 
Bxy ~ Bt sign{z) sin 6 , 
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/iT = 1.5 kpc is the height of the maximum above the 
plane and wt = 0.3 kpc is its lorentzian width. In con- 
trast to HJI, we choose 



ST,max(r) = l.biiG Q{Ro -r) + e(r - i?o)e— 

' (9) 

with 8 denoting the Heaviside step function, so that the 
halo contribution becomes negligible for r ^ Rq, because 
there is no evidence for such a field outside the solar 
circle Finally, a dipole field is added as e.g. in 0, 
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(10) 



4 



where R = ^r'^ + z'^ = a/x^ + + z^, cos (p = z/R and 
/iG is the magnetic moment of the Galactic dipole with 
/iG = 123/.jGkpc'^ in order to reproduce — +0.2^G 
near the Solar system . To avoid a singularity in the 
center, we set = — 100/^G inside a sphere of 500 pc 
radius centered at the GC. Note that in [s^l values as 
large as 1 mG were used for the hard core of the dipo- 
lar field. However, data from low frequency non-thermal 
radio emissions of electrons [s^l favor a value of about 
10 juG down to a 10 pc scale, and put the safe bound of 
100 ^G which we actually use. 

Finally, we warn the reader that these models intend to 
provide only a rough approximation to the true structure 
of the GMF. Moreover, the GMF models are not self- 
consistent because the condition V • B = is not fulfilled 
by any of the disk fields discussed: using the ansatz of 
Eq. Q for S(r, 6*), V • B can only be satisfied by 
h{r) = const, or a non-vanishing z-component of the disk 
field. 



III. GMF AND UHECR PROPAGATION 

A generalized version of the Liouville theorem was 
shown to be valid for CRs propagating in magnetic fields 
soon after the discovery of the geomagnetic effect [36j . 
For UHECRs, this theorem has been numerically tested 
e.g. in [33, where particles were injected isotropically 
from randomly distributed sources at different Galac- 
tocentric distances. Even after the propagation in the 
GMF, the sky on Earth appears isotropic (see their Fig. 6, 
left panel), while the effective exposure of an experiment 
to the extragalactic sky is strongly modified by the GMF 
(see their Fig. 6, right panel). Indeed, whenever the rigid- 
ity cutoff can be neglected for high energy cosmic rays, for 
an isotropic fiux outside the Galaxy the GMF introduces 
no anisotropy. At sufficiently low energy, however, the 
GMF may introduce blind regions on the external sky, 
which translate into observed anisotropics for an Earth- 
based observer. This effect is easy to estimate in the 
simple case in which our Galaxy has a large scale dipo- 
lar field as the one introduced in the PS model. The 
St0rmer theory (sec e.g. [|[) can be applied to determine 
the rigidity cutoff TZs below which no particle can reach 
the Earth. Since the Earth is at zero Galacto-magnetic 
latitude, we obtain 
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Here, TZs depends on the arrival direction of the cosmic 
ray via the function e, that we do not need to specify here. 
Assuming that the tiny vertical component detected at 
the solar system of 0.2 jiG is due to a dipole field, we get 
fiQ — 120 /xG kpc'^ and then TZs would vary in the range 
10^^ V-10^^ V. For such a GMF, large-scale anisotropics 
should be seen around E ~ ZeTZs, if an extragalactic 
component dominates at this energy. Models that invoke 
a dominating extragalactic proton component already at 



E ~ 4 X 10^^ eV (see e.g. jSg) or extragalactic iron nu- 
clei at £^ < 10^^ eV might be in conflict with the observed 
isotropy of the CR flux. Unfortunately, the previous con- 
clusion is not very robust against changes in the GMF 
model. Indeed, the geometry of the GMF is more com- 
plicated than a simple dipole. A naive estimates of the 
rigidity cutoff TZs might be obtained from the size of the 
Larmor radius in the Galaxy, 
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(12) 



For B ~ few fiG and a Galactic magnetic disk of thick- 
ness 0(100) pc, Eq. nil) implies 7^ < few x lO^'^ V. 
However, the O(IOO) pc disk thickness is comparable to 
the largest scale of turbulent fields. Thus, quasi-diffusive 
processes might work in an opposite direction, partially 
erasing anisotropics due to blind region effects in an ex- 
tragalactic dominated fiux. One should also keep in mind 
that in the same energy range a sizeable anisotropy is 
expected in Galactic models of the origin of CRs, al- 
though usually with a different large scale pattern [3^ . 
Given present knowledge of the GMF, we conclude that 
the signature of blind regions for an extragalactic fiux at 
TZ $ few X 10^^ V is less robust and unambiguous than 
naively estimated. 

The previous discussion referred to the case of an in- 
coming isotropic flux. However, if the AGASA small- 
scale clusters are not just a statistical fluctuation, the 
UHECR flux is, at least on small-scales, anisotropic. In 
this case, the CR flux can be magnified or de-magnified 
by magnetic lensing phenomena, and the application of 
the Liouville theorem is non-trivial |40j . The magnifica- 
tion effects of the GMF change the experimental expo- 
sure and a well-defined procedure is needed to assess the 
significance of any detected anisotropy. One can test the 
significance of observed anisotropics by comparing the 
values of the statistical estimator based on the Nd data 
with a large number Af of simulated iV^-points samples 
of the null-hypothesis. For each set, one should consider 
the propagation in the GMF, convolve with the exper- 
imental exposure, and finally reject the null-hypothesis 
with a given significance. In the following, we use for 
practical reasons the backtracking technique Since 
we deal with ultra- relativistic particles, the equations of 
the motion can be written in the form 



dv V X B 

It ~ TZ 



(13) 



where v is a vector of modulus practically equal to c = 1. 
The integration is stopped when the particle reaches a 
distance of 50 kpc from the GC. Note that the energy 
losses of UHECRs on Galactic scales (^ 10 kpc) are neg- 
ligible, provided the trajectory is not very far from a 
rectilinear one. In Fig. we show the map thus ob- 
tained, for the three models considered and a rigidity of 
4 X 10^^ V. Note that the deficction (5 of a particle of 
rigidity E/Ze in a field of strength B coherent over the 
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FIG. 2: Deflection maps for the TT (top), HMR (middle), 
and PS (bottom) models of the GMF, for a rigidity of 4x 
10^^ V. The deflection scale is in degree, and the maps refer 
to the direction as observed at the Earth. All maps use a 
Hammer-Aitoff projection of Galactic coordinates. 



scale L is approximately given by 
S 



„ Z B L 
0.53°- - — 

ii20 fJ-(j kpC 



(14) 



To obtain an estimate for the average deflection of 
CRs in different sky regions, labeled as in Table ^ we 
have followed backwards 50000 randomly chosen CRs of 
rigidity TZ ~ 10^" V, for which the hypothesis of a quasi- 
rectilinear trajectory is well fulfilled. In Tablc^I wc show 



the average value and the square root of the variance of 
7?.2o(5 (i.e., in units of lO^*' V) for the three models of 
the regular GMF discussed in the Section calculated 
separately for eight different sky regions. The quantity 
TZ6 depends only on the GMF model and scales almost 
linearly with the field strength B. The largest difference 
between the three GMF models occurs in the region Al: 
in the PS model, the only one with a dipole field, huge 
deflections arise close to the GC. In the regions Bh, Ch, 
and Dh the stronger halo fields of the TT and especially 
the HMR model cause larger defiections than in the PS 
model. In the 1-regions, apart for Al, the defiections of 
the three models are all of the order 1°~2°, and compa- 
rable to each other within 1 <t. Since the CR in these 
directions mainly travel through the disk, in order to es- 
cape the Galaxy they have to cross the regions where 
the field geometry and intensity is better known, and a 
better agreement among the models exists. 

If one excludes the central regions of the Galaxy, the 
average deflection is (5 ~ 2°/7?.2o, and the differences for 
the magnitude of the deflections are of the order of 50% 
among the models. Thus only for the highest energy 
events and proton primaries the role of the GMF is neg- 
ligible compared to the angular resolution Scxp of CR 
experiments. The latter is as good as (5oxp — 0.6° for the 
HiRcs experiment Q and for Auger hybrid events [i^ . 
For lower TZ, correcting for defiections in the GMF would 
becomes crucial to exploit fully the angular resolution of 
UHECR experiments. Note also that, even in the ideal 
case of perfectly known GMF, a reconstruction of the 
original arrival directions would require a relatively good 
energy resolution: an uncertainty of, say, 30% in the en- 
ergy scale around 5 x 10^^ eV would lead to errors ^ 1° 
in the reconstructed position of proton primaries in most 
of the sky. 

Apart for deflections, (de-) focusing effects of the GMF 
effectively modify the exposure to the extragalactic sky. 
This modification can be calculated by back-tracking a 
large number of events and looking at the obtained map 
of event numbers per solid angle outside the Galaxy. For 
the purpose of illustration, we show in Fig. |31 some rel- 
ative exposure maps obtained for fixed rigidity for the 
three chosen GMF models. They were obtained with the 
technique described in |33j . and essentially represent the 
ratio LOB{l,b) = dQoail' ,b')/dfli^{l,b), where dfl^ is an 
infinitesimal small cone at Earth (around the direction 

b) transported along the trajectory of a charged par- 
ticle to the border of the Galaxy d^oc (around the new 
position l'{l,b),b'{l,b)). If ujB{l,b) deviates significantly 
from one, the corrected exposure has to be used in (auto- 
) correlation studies. Note how this effect is present in 
all the models at least for cosmic rays observed at the 
Earth along the Galactic plane. 

A remark on the role of the turbulent GMF is in or- 
der. A comparison of Eq. I^J and Eq. H14() shows that 
defiections in the turbulent GMF are sub-leading. How- 
ever, this does not ensure automatically that magnetic 
lensing by the turbulent fields is irrelevant at high cner- 
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baVc 


315 < Ig <45 


45<Ig < 135 


135 <Ig < 225 


225 <Ig < 315 


[bal >30 


Ah 


Bh 


Ch 


Dh 


I&gI <30 


Al 


Bl 


CI 


Dl 



TABLE I: The labels used for the eight different regions of the sky referred to in the text. Angles are in degrees. 



region 


7^2o s (PS) 


7^2o S (TT) 


7^2o S (HMR) 


Ah 


1.8° ± 1.0° 


0.9° ±0.5° 


2.1° ±0.8° 


Bh 


1.3° ± 0.4° 


1.3° ± 0.6° 


2.2° ±0.8° 


Ch 


0.9° ± 0.4° 


2.0° ±0.3° 


2.7° ±0.5° 


Dh 


0.5° ± 0.2° 


1.1° ± 0.6° 


2.1° ±0.9° 


Al 


14° ±21° 


1.9° ±0.4° 


2.2° ±0.7° 


Bl 


2.0° ± 0.9° 


1.7° ± 0.5° 


1.2° ±0.4° 


CI 


1.7° ± 1.1° 


1.9° ± 0.5° 


1.8° ±0.3° 


Dl 


1.6° ± 1.0° 


1.6° ± 0.5° 


2.3° ±0.6° 



TABLE IL The rigidity times average deflections 72.20 3 in the eight different regions labeled in Table |I| 



gies, because lensing depends on the gradient of the field 
and the critical energy for amplification is proportional to 
^-1/2^ The detailed analysis of showed indeed that 
small-scale turbulence can produce relatively large mag- 
nification effects, and argued that it may even be respon- 
sible for (some of) the multiplets seen by AGASA above 
4 X 10^^ eV. However, we expect that the random GMF 
introduces a small scale (and strongly energy dependent) 
correction of the sky map on the top of the magnification 
effects of the regular GMF. Then, we argue that neglect- 
ing the turbulent GMF in our analysis is justified both 
for intrinsic reasons and experimental limitations. First, 
possible lensing effects by the random GMF are weak- 
ened by the presence of a regular field component |l5j |: 
Already without regular field, the magnification peaks 
are quite narrow in energy space, AE/E ^ 20%, and 
thus their width AE/E is comparable to the energy res- 
olution of CR experiments. The presence of the regu- 
lar GMF narrows these magnification peaks further by a 
factor of several. Due to the limited energy resolutions 
of current experiments, the turbulent lensing signature 
would be smeared out. Second, the experimental angu- 
lar resolution introduces an additional averaging effect. 
Thus it is seems unlikely that lensing effects of the tur- 
bulent GMF lead at present to distinctive effects taking 
into account the current experimental limitations. Fi- 
nally, note that from a phenomenological point of view 
our analysis is the same both if the multiplets are due to 
intrinsically strong UHECRs sources or due to turbulent 
lensing. 



IV. AGASA DATA SAMPLE 

In order to make the general considerations of the pre- 
vious section more concrete, we will discuss here some 



applications to the AGASA data. The AGASA exper- 
iment has published the arrival directions of their data 
until May 2000 with zenith angle < 45° and energy above 
4 X IQis eV HI. This data set consists oi N = 57 CRs 
and contains a clustered component with four pairs and 
one triplet within 2.5° that has been interpreted as first 
signature of point sources of UHECRs. The reconstruc- 
tion of the original CR arrival directions and the estimate 
of their errors is obviously an important first step in the 
identification of astrophysical CR sources. 

In Fig. 21 we show the measured directions of all CRs 
in the AGASA data with E > 10^° eV together with 
their reconstructed arrival direction at the border of the 
Galaxy for the case of Carbon primaries. In the south- 
ern Galactic hemisphere, the TT model often produces 
opposite deflections with respect to the HMR and PS 
models, because of its antisymmetric field configuration. 
A longitudinal shift is often appreciable in the PS model, 
as consequence of the dipolar component we added. The 
larger deflection found in the TT and HMR models for 
high latitudes is explained by the stronger regular halo 
field these models use. The best chances for source iden- 
tification arise obviously by looking at directions opposite 
to the GC. On the other hand, observations towards the 
GC have a certain importance to use the UHECRs as 
diagnostic tool for GMF, see e.g. [i^ . 

Before turning to the statistical analysis, we briefly 
recall the estimators we will use in the following. The 
autocorrelation function wi is defined as 

i=l j=l 

where lij is the angular distance between the two cosmic 
rays i and j, £ the chosen bin size, O the step function, 
and Nd is the number of CRs considered. Analogously, 



7 



0.1 1 10 




FIG. 3: The extragalactic "exposure" maps for the TT (top), 
HMR (middle), and PS (bottom) model for a rigidity of 4x 
10^^ V. AU maps use a Hammer- AitofT projection of galactic 
coordinates. 



one can define the correlation function Wg as 

i—1 a—1 

where iia is the angular distance between the CR i and 
the candidate source a and Ng is the number of source 
objects considered. 
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FIG. 4: The AGASA data set of CRs with estimated energy 
exceeding 10^° eV together with their reconstructed positions 
assuming Carbon primaries and the three GMF models dis- 
cussed in the text. The line bounds the region visible to 
AGASA for zenith angles < 45°; note that for a better visi- 
bility, this map is centered around / = 180°, since AGASA is 
blind to the Galactic Center. 



A. Autocorrelation analysis 

Let us discuss now how the small-scale clustering ob- 
served by the AGASA experiment is modified by the 
GMF. Note from Fig. 13 that even for protons the (de- 
) magnification of the exposure is already significant at 
energies 4x 10^^ eV. Neglecting the influence of the GMF 
(i.e., assuming neutral primaries), one generates a large 
number of Monte Carlo sets of CRs, each consisting of 
Nd CRs distributed according to the geometrical expo- 
sure Woxp of AGASA. The fraction of MC sets that has 
a value of the first bin of the autocorrelation function 
wi larger or equal to the observed one, wl, is called 
the chance probability P of the signal. In the follow- 
ing, we shall consider a bin- width of 2.5°, suited for 
AGASA angular resolution. For a nonzero GMF, one 
uses the back-tracking method: the observed arrival di- 
rections on Earth are back-tracked following a particle 
with the opposite charge to the boundary of the GMF. 
Then the value w of the autocorrelation function is cal- 
culated. Since also the exposure is changed by the GMF, 
the CRs of the Monte Carlo sets have to be generated 
now using as exposure uJtoti'R-, l,b) = Wexp(^, &)^^b(7?., h b) 
(the energy spectrum is sampled according to the spec- 
trum observed by AGASA, d7V/d£: oc E''^-'^). This is 
automatically fulfilled if one back-tracks uniformly dis- 
tributed Monte Carlo sets in the same GMF^. The re- 



For very large statistics, however, it is numerically more conve- 
nient to explicitly calculate u)B{T^,l,b) and generate the Monte 
Carlo sets accordingly. 
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suiting chance probability is called P in Table Hill The 
error due to the statistical fluctuation in the Monte Carlo 
sets is on the last digit reported. For illustration, we show 
also the chance probability Pq calculated using only the 
experimental exposure (or lub = 1) that overestimates 
that clusters come indeed from the same source. 

Correcting for the GMF reduces for all three GMF 
models the value of wi. While however for the TT and 
HMR models two doublets above 5 x 10^^ eV disappear, 
the PS model loses one low-energy doublet. Thus, either 
some of the pairs are created by the focusing effect of the 
GMF, or the GMF and especially its halo component is 
not well enough reproduced by the models. In the lat- 
ter case, "true pairs" are destroyed by the the incorrect 
reconstruction of their trajectories in the GMF. Refer- 
ence ^3 discussed in detail the effect of the GMF on 
the AGASA triplet and found that current GMF mod- 
els defocus it. If the clustering is physical, this could 
be explained by a wrong modeling of the GMF in that 
high-latitude region. Alternatively, our assumptions of 
negligible deflections in the extragalactic magnetic field 
and of protons as primaries could be wrong. Note that 
the effect of the GMF induced exposure to the extra- 
galactic sky is not in general negligible. The fact that 
Pq is only somewhat smaller than P hints that only a 
small fraction of clusters might be caused by magnetic 
lensing (in the regular field). For the AGASA data set 
this is expected, because the field of view of this experi- 
ment peaks away from the inner regions of the Galactic 
plane, and the GC in particular. Finally, we note that 
the energy threshold for which the chance of clustering 
is minimal decreases for the TT model. This change is 
however rather small and a larger data set is needed for 
any definite conclusion. 

In Table IIVI the same analysis is performed for the 
TT model only, but assuming also Z = +2, —1. In no 
case is the improvement with respect to the Z = case 
appreciable. 



B. Correlations with BL Lacs 

Tinyakov and Tkachev examined in p^ ] if correlations 
of UHECRs arrival directions with BL Lacs improve af- 
ter correcting for deflections in the GMF (see also ^3). 
The signiflcance of the correlation found is still debated, 
and we just choose this example as an illustration how 
correlation of UHECR arrival directions with sources can 
be used to test the GMF model and the primary charge. 
We use from the BL Lac catalogue ^3 a-H confirmed BL 
Lacs with magnitude smaller than 18 (187 objects in the 
whole sky). 

In Table EI we show the chance probability to observe 
a stronger correlation taking into account the three dif- 
ferent GMF models and assuming proton primaries. An 
improvement of the correlation signal is found only for 
the TT model, while for the two other models the corre- 
lation becomes weaker. In Table fVll the same analysis is 



performed for Z = +2 and —1 using the TT model. An 
improvement with respect to the Z = case is only found 
for protons, i.e. for Z = +1. Obviously, the two cases 
Z = +2 and -1, i.e. a CR flux consisting only of Helium 
nuclei or anti-protons, is not expected theoretically, and 
is used mainly to illustrate the potential of this method 
to distinguish between different charges. Although it is 
difficult to draw strong conclusions given present uncer- 
tainties, this example clearly shows how UHECRs ob- 
servations might be used to constrain the GMF models 
and/or to determine the charge of the CR primaries. 



V. CONCLUSION AND PERSPECTIVES 

In this work we have discussed in detail the effect of the 
regular component of the Galactic magnetic field on the 
propagation of UHECR, using three GMF models dis- 
cussed previously in the literature. Both in small-scale 
clustering and correlation studies, the GMF might be 
used as a sort of natural spectrograph for UHECR, thus 
helping in identifying sources, restricting the GMF mod- 
els as well as the chemical composition of the primaries. 
Notice that the latter point is an important prerequisite 
to use UHECR data to study strong interaction at energy 
scales otherwise inaccessible to laboratory experiments. 

As a consequence of the existence of blind regions, 
we have argued that in some GMF models the observed 
isotropy of cosmic ray fiux e.g. around a few 10^^ eV 
might disfavor a transition to an extragalactic flux of 
protons at too low energies, say below 10^^'^ eV. Unfor- 
tunately the effect is model-dependent, and a better char- 
acterization of the GMF is needed to draw more robust 
conclusions. 

At higher energies, the cosmic rays should enter the 
ballistic regime. We have provided in tabular and picto- 
rial form an estimate of the typical deflection suffered by 
UHECRs in the small deflection limit and of its variabil- 
ity from model to model. These results imply that, if the 
angular resolution of current experiments has to be fully 
exploited, deflections in the GMF cannot be neglected 
even for E = 10^° eV protons, especially for trajectories 
along the Galactic plane or crossing the GC region. Since 
the magnitude as well as the direction of the deflections 
are very model-dependent, it is difficult to correct for 
deflections with the present knowledge about the GMF. 

We have also emphasized that, to the purpose of sta- 
tistical analyses like auto-correlation/cross-correlation 
studies, the GMF can effectively act in distorting the 
exposure of the experiment to the extragalactic sky, and 
we provided some maps of this "exposure-modification" 
effect. As a consequence, to estimate the chance proba- 
bility e.g. of small-scale clustering, one should take into 
account the role of the GMF. We showed that this effect 
is already appreciable in the data published by AGASA, 
although its fleld of view do not include the central 
regions of the Galaxy. Especially for experiments in the 
southern hemisphere like Auger, one might wonder if 
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Model 


no GMF 


TT 


HMR 


PS 


£n.in/10^^ eV 


Nd 


Wi P[%] 


w, P[%] Po[%] 


Wi P[%] Po[%] 


P[%] Pq[%] 


5.0 
4.5 
4.0 


32 
43 
57 


4 0.22 

6 0.05 

7 0.18 


2 8.9 8.5 
4 1.5 1.4 
6 0.8 0.7 


2 9.5 8.5 

4 1.8 1.4 

5 3.3 2.4 


4 0.27 0.22 
6 0.10 0.05 
6 1.09 0.72 



TABLE III: Number A'' of CRs with energy E > -Emin and zenith angle 9 < 45°; the values of the first bin of the autocorrelation 
function wi, and the chance probability P{wi > w*) from an isotropic test distribution are shown for the two cases with (P) 
and without (Po) correction of the exposure, respectively. Proton primaries are assumed. 



Model 




Z 


=+1 


Z 


=+2 


Z=-l 


£;mi„/10^^ eV 


Ad 


Wi P[%] 


Wi 


P[%] 


Wi 


P[%] 


P[%] 


5.0 


32 


4 0.22 


2 


8.9 


1 


41 


4 0.27 


4.5 


43 


6 0.05 


4 


1.5 


4 


1.8 


6 0.10 


4.0 


57 


7 0.18 


6 


0.8 


5 


3.5 


6 1.08 



TABLE IV: As in Tabic ETTI but for diflterent charges of the primary (TT model). 



Model 


no 


GMF 


TT 


HMR 


PS 


Sn.in/lO^*' eV 


Ad 




P[%] 




P[%] 


Po[%] 




P[%] Po[%] 


■Wbl 


P[%] 


Po[%] 


5.0 


32 


8 


8.8 


8 


9.7 


8.8 


4 


65 64 


7 


17 


17 


4.5 


43 


11 


4.7 


14 


0.6 


0.5 


7 


49 39 


7 


39 


39 


4.0 


57 


13 


6.6 


20 


0.05 


0.05 


11 


20 18 


7 


67 


67 



TABLE V: As in Table UTTl but for the correlation function and the BL Lac catalogue discussed in the text. 



Model 


Z=0 


Z-- 


=+1 


Z-- 


=+2 


Z 


=-1 


B^^in/lO^^ eV 


Ad 


Wbi P[%] 




P[%] 


■Whl 


P[%] 




P[%] 


5.0 


32 


8 8.8 


8 


9.7 


7 


20 


6 


27 


4.5 


43 


11 4.7 


14 


0.6 


12 


3.4 


7 


36 


4.0 


57 


13 6.6 


20 


0.05 


13 


9.9 


8 


50 



TABLE VI: As in Table El but for different charges of the primary (TT model). 



excluding some part of the data from (auto-) correlation 
studies might lead to more robust analyses, at least as 
long as no reliable model for the GMF is established. 
We stress that the required cuts would be quite drastic. 
For instance, fixing i^min = 4 x 10^^ eV and consid- 
ering only sky regions where jt^s — 1| < 0.2 would exclude 

TT: -5° < 5g < 5° for all la and -60° < ba < 60° 
for -94° <Ig< 75°, 

PS: -25° < be < 22° for all h and -38° < be < 40° 
for -33° <Ig< 35°, 

HMR: -11° < ba < 8° for aU Iq and all be for 
-90° <Ig < 90°. 

Note that in the HMR model one would cut more 
than half of the sky. A more reasonable prescription is 
to evaluate the robustness of the significance of autocor- 
relation and correlation studies by taking into account 



several models of GMF and several primary charges, 
as we have done in this paper for the public available 
AGASA catalogue. Although the present statistics does 
not allow to draw strong conclusions, we have not found 
any signal of improvement after the correction for GMF. 
This could point to an insufficient knowledge of the field 
or to a significantly heterogeneous chemical composition 
of the primaries. Finally, the AGASA signal might only 
be a chance fluctuation. 

Independently on the outcome, performing statistical 
analyses taking into account several models of GMF and 
several primary charges is an useful exercise. In the most 
pessimistic case, it allows one to quantify in an approx- 
imate way the contribution of the GMF to the overall 
uncertainty in the chance probability of a candidate sig- 
nal. On the other hand, a strong improvement in the sig- 
nificance of a statistical estimator might favor a certain 
GMF model and/or primary charge assignment. For ex- 
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ample, by repeating the study of |2J| we have found that 
the significant correlation of BL Lacs with UHECRs is 
strongly dependent on the GMF and primary adopted, 
and is present only in the TT model of the GMF for 
Z = \. Although this evidence needs confirmation with 
a larger data set of UHECRs, it may be the start of the 
era of UHECR astronomy. 



comments on the manuscript. M.K. was partially sup- 
ported by the Deutsche Forschungsgemeinschaft within 
the Emmy-Noether program. PS acknowledges the 
support by the Deutsche Forschungsgemeinschaft under 
grant SFB 375 and by the European Network of The- 
oretical Astroparticle Physics ILIAS/N6 under contract 
number RII3-CT-2004-506222 



Acknowledgments 



We are grateful to Dmitry Semikoz for collaboration 
during the initial phase of this work and for useful 



[1] For recent reviews see e.g. M. Nagano and A. A. Wat- 
son, Rev. Mod. Phys. 72, 689 (2000); J. W. Cronin, 
^stro-ph/0402487 M. Kachelriefi, Comptes Rendus 
Physique 5, 441 (2004) hep-ph/0406174 . 

[2] A. A. W atson, Nucl. Phys. Proc. Suppl. 136, 290 (2004) 
jastro-ph/0408110 . 

[3] roT an early review see Th. H. Johnson, Rev. Mod. Phys. 
10, 193 (1938). 

[4] H. Takami, H. Yoshiguchi and K. Sato, Astrophys. J. 
639, 803 (2006) astro- ph/0506203 . 

[5] Y. Uc hihori et al, Astropart. Phys. 13, 151 (2000) 
|astro- ph/9908193 . M. Takeda et al, Astophys. J. 522 
225~'astro-ph/9902239 ; M. Takeda et al, Proc. 27th 
ICRC, Hamburg, 2001. 

[6] P. G. Tinyakov and I. I. Tkachev, JETP Lett. 74, 
1 (2001) [Pisma Zh. Eksp. Teor. Fiz. 74, 3 (2001)] 
^stro-ph/0102101j; C. B. Finley and S. Westerhoff, As- 
tropart. Phys. 21, 359 (2004) astro-ph/0309159 . 

[7] H. Yoshiguchi, S. Nagata ki and K. Sato, Astrophys. J. 
614, 43 (2004) "astro-ph/0404411 . 

[8] R. U. Abbasi et al |The High Resolution Fly's Eye 
Collaboration (HIRES)], Astrophys. J. 610 (2004) L73 
rastro-ph/0404137 . 

[9] M. Kachelriefi and D. Semikoz, Astropart. Phys. 23, 486 
(2005) astro- ph/0405258 . 
[10] B. Revenu et al [The Pierre Auger collaboration], Proc. 

29th ICRC, Pune, 2005. 
[11] K. Dolag, D. Grasso, V. Springel and I. Tkachev, 
JETP Lett. 79, 583 (2004) [Pisma Zh. Eksp. Teor. 
Fiz. 79, 719 (2004)] astro-ph/0310902 ; JCAP 0501, 
009 (2005) astro- ph/0410419 ; G. Sigl, F. Miniati 
and T. Enfilin, Phys. Rev. D 70, 043007 (2004) 
astro-ph/0401084 ; Nucl. Phys. Proc. Suppl. 136, 224 

(2004) astro- ph/0409098 . 

[12] W. A. Hiltner, Astrophys. J. 109, 471 (1949). 

[13] E. Zwiebel and C. Heiles, Nature 385, 131 (1997). 

[14] R. Stepanov, P. Frick, A. Shukurov and D. Sokoloff, As- 

tron. Astrophys. 391, 361 (2002) astro-ph/0112507 . 
[15] D. Harari, S. MoUerach, E. Roulct and F. Sanchez, JEEP 

0203, 045 (2002) astro- ph/02 02362 . 
[16] P. G. Tinyakov and I. I. Tkachev, Astropart. Phys.24, 32 

(2005) astro- ph/041 1669 . 

[17] D. Sokolov, A. Shukurov and M. Krause, Astron. Astro- 
phys. 264, 396 (1992). 
[18] P. Frick, R. Stepanov, A. Shukurov and D. Sokoloff, 



Mon. Not. Roy Astron. Soc. 325, 649 (2001) 
astro-ph/0012459 . 
[19] R. Beck, in "The Astrophysics of Galactic Cosmic Rays", 
eds. R. Diehl et al. Space Science Reviews, Kluwer 
astro-ph/0012402 . 
[20] R. Beck, A. Brandenburg, D. Moss, A. Shukurov, 
D. Sokoloff, Ann. Rev. Astron. Astrophys. 34, 155 
(1996). 

[21] R. J. Wainscoat et al, Astrophys. J. Suppl. 83, 111 
(1992). 

[22] J. Han, in "Astrophysical Polarized Backgrounds", 

American Institute of Physics, edited by S. Cecchini et 

al, p. 96 (2002) astro-ph/0110319 . 
[23] T. Stanev, Astrophys. J. 479, 290 (1997). 
[24] P. G. Tinyakov and I. I. Tkachev, Astropart. Phys. 18, 

165 (2002) astro-ph/0111305 . 
[25] R. J. Rand and A. G. Lync, Mon. Not. Roy. Astron. Soc. 

268, 497 (1994). 
[26] K. Beuermann, G. Kanbach, E. M. Berkhuijsen, Astron. 

Astrophys. 153, 17 (1985). 
[27] J. L. Han, R. N. Manchester, E. M. Berkhuijsen, R. Beck, 

Astron. Astrophys. 322, 98 (1997). 
[28] J. L. Han, R. N. Manchester, G. J. Qiao, Mon. Not. Roy. 

Astron. Soc. 306, 371 (1999). 
[29] J. L. Han and G. J. Qiao, Astron. Astrophys. 288, 759 

(1994). 

[30] J. L. Han, R. N. Manchester, G. J. Qiao and A. G. Lyne, 

in "R adio Pulsars", APS Conf. Ser., editted by M. Bailes 

et al |astro-ph/0211197| . 
[31] S. Roy, Bull. Astron" Soc. India 32, 205 (2004) 

astro- ph/041 1687 . 
[32] T. N. LaRosa et al, Astrophys. J. 626, L23 (2005) 

astro- ph/0505244 . 
[33] D. Harari, S. Mollcrach and E. Roulet, JHEP 9908, 022 

(1999) astro-ph/9906309 . 
[34] M. Prouza and R. Smida, Astron. Astrophys. 410, 1 

(2003) astro-ph/0307165 . 
[35] H. Yoshiguchi, S. Nagataki and K . Sato, Astrophys. J. 

596, 1044 (2003) astro-ph/0307038 . 
[36] G. Lemaitre and M. S. Vallarta, Phys. Rev. 43, 87 (1933); 

W. F. G. Swann, Phys. Rev. 44, 224 (1933). 
[37] J. Alvarez-Muniz, R. Engel and T. Stanev, Astrophys. J. 

572, 185 (2002) astro-ph/0112227 . 
[38] M. Ahlers et al, Phys. Rev. D 72, 023001 (2005) 
astro-ph/0503229, . 



11 



[39] J. Candia, S. Mollerach and E. Roulet, JCAP 0305, 003 

(2003) astro- ph/0302082 . 
[40] D. Harari, S. Mollerach and E. Roulet, JHEP 0002, 035 

(2000) astro-pli/0001084 . 
[41] S. Karakula, J. .L. Osborne, E. Roberts and W. Tkaczyk, 

'Proc. 12th ICRC (Hobart) 1, 310 (1971). 
[42] C. Bonifazi et al. [The Pierre Auger collaboration], Proc. 

29th ICRC, Pune 2005. 
[43] N. H ayashida et al, Astron. J. 120, 2190 (2000) 

|astro -ph/0008102 . 
[44] D. Harari, S. Mollerach and E. Roulet, JHEP 

0207, 006 (2002) astro-ph/0205484 G. Medina-Tanco, 

M. Teshima, and M. Takeda, Proc. 28th ICRC, Tsukuba 

2003; J. Alvarez-Muniz and T. Stanev,|astro-ph/0507273| 



[45] S. V. Troitsky, astro-ph/0505262 

[46] For a discussion of possible correlations of BL Lacs with 
the HiRes data see D. S. Gorbunov et al, JETP Lett. 80, 
145 (2004) [Pisma Zh. Eksp. Teor. Fiz. 80, 167 (2004)] 
astro- ph/0406654 ; R. U. Abbasi et al [HiRes Collabora- 
tion], Astrophys. J. 636, 680 (2006) astro-ph/0507120|j. 
The correlation signal between UHECRs and BL Lacs in 
futute data sets is estimated in D. S. Gorbunov et al, 
JCAP 0601, 025 (2006) astro-ph/0508329 . 

[47] M.-P. Veron-Cetty and P. Veron, Astron. Astrophys. 
412, 399 (2003). The catalogue is available on-line at 
http: //www. obs-hp .f r/www/ catalogues/veron2_ll/veron2_ll .htm 



